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Bundles of many fibers, with statistically distributed thresholds for breakdown of individual fibers 
and where the load carried by a bursting fiber is equally distributed among the surviving members, 
are considered. During the breakdown process, avalanches consisting of simultaneous rupture of 
several fibers occur, with a distribution D(A) of the magnitude A of such avalanches. We show that 
there is, for certain threshold distributions, a crossover behavior of D(A) between two power laws 
D(A) oc A -4 , with £ = 3/2 or £ = 5/2. The latter is known to be the generic behavior, and we give 
the condition for which the D(A) oc A~ 3 ^ 2 behavior is seen. This crossover is a signal of imminent 
catastrophic failure in the fiber bundle. We find the same crossover behavior in the fuse model. 

PACS numbers: 



A fundamental question in strength considerations of 
materials is when does it fail. Are there signals that can 
warn of imminent failure? This is of uttermost impor- 
tance in e.g. the diamond mining industry where sudden 
failure of the mine can be extremely costly in terms of 
lives. These mines are under continuous acoustic surveil- 
lance, but at present there are no tell-tale acoustic signa- 
ture of imminent catastrophic failure. The same type of 
question is of course also central to earth quake predic- 
tion. We will in this letter study signatures of imminent 
failure in the context of the fiber bundle model We 
find that if a histogram of the number of fibers failing 
simultaneously is recorded over an interval which starts 
sometime during the failure process, it follows a power 
law with an exponent that crosses over to a very dif- 
ferent value if the start of the interval is close enough 
to the point at which the fiber bundle fails catastrophi- 
cally. This is a clear signature of imminent failure. We 
also study the fuse model in this context and find that 
it behaves qualitatively in the same manner as the fiber 
bundle model. 

When a weak element in a loaded material fails, the in- 
creased stress on the remaining elements may cause fur- 
ther failures, and thereby give a burst avalanche in which 
A elements fail simultaneously. With further increase 
in the load new avalanches occur. A bundle of many 
fibers with stochastically distributed fiber strengths, and 
clamped at both ends, is a much studied model for such 
avalanches. In its classical version a ruptured fiber 
carries no load and the increased stresses caused by a 
failed element are shared equally by all surviving fibers. 

A main result 0iQ for this model is that under mild re- 
strictions on the fiber strength distribution the expected 
number D(A) of burst avalanches of size A is governed 
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FIG. 1: The distribution of bursts for the strength distribu- 
tion @ with xo = and xq = 0.45x m . The figure is based on 
10000 samples with N = 50000 fibers. 



by a universal power law 

£>(A) cx A - * 



(1) 



for large A, with £ = 5/2. 

When the load on the bundle is increased beyond a crit- 
ical threshold x c , the whole bundle ruptures 0,0]. We 
are interested in the final stages of the breakdown pro- 
cess, when the N surviving fibers has a distribution of 
fiber strengths in the interval {xq, x m ), with xq slightly 
less than x c ). Assume first that the fiber strengths in 
which the maximal loads x n that the fibers n = 1, 2, . . . N 
are able to carry are picked independently with a proba- 
bility density 

p(x) — Prob(a; < x n < x + dx) j dx 
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(x m — x ) 1 for Xq < x < x r 



(2) 



When xq is greater than the critical value x c = x m /2 the 
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FIG. 2: The distribution of bursts for the strength distribu- 
tion (Eqn.|^l with xo — 0A0x m - The figure is based on 10000 
samples with N = 50000 fibers. The straight lines represent 
two different power laws and the arrow locates the crossover 
point A c ~ 12.5. 



whole bundle fails fatally when the load per fiber exceeds 
xq. This can be seen as follows: The force F(x) that the 
bundle is able to withstand when all fibers with strengths 
less than x have ruptured, is proportional to the number 
of surviving fibers times the strength, 



F(x) (xx J p{x) dx 

for xq < x < x m . 
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When xo > x c = x m /2, this is a decreasing function of x. 
The consequence in this case is that after the first fiber 
ruptures at x = Xq the whole bundle will fail completely 
at once 5] . A threshold distribution in which the weakest 
fiber has the critical value we call a critical threshold 
distribution. We want to study the burst distribution 
close to this critical situation. In Fig. ^ we show results 
for -D(A) in a simulation experiment on fiber bundles 
with the strength distribution with xo = 0A5x m — 
0.9x c . For comparison, results with xq = are shown. 
In both cases -D(A) shows a power law decay, apparently 
with an exponent £ = 5/2 for the xq = case, and an 
exponent 3/2 for xq = 0A5x m (bundle with no weak 
fibers) 0. In this note we explain the latter result as a 
crossover phenomenon. 

For a bundle of many fibers the expected number of 
bursts of magnitude A is given by 
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-J Xc p(x) [l-xp{x)/Q{x)] [xp(x)/Q(x)] A - 1 
x exp [—Axp(x)/Q(x)] dx , 



(4) 



where Q(x) = p(y) dy is the fraction of fibers with 
strength exceeding x. For the distribution J2J this yields 
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Introducing the parameter 

X c ~ X 

e = 

and a new integration variable 

Xm 

- 2x 



we obtain 

D(A) _ 2A A ~ 1 e- A e 2 
N ~ A!(l + 2e) 

X JO (l-ez)(2-tz) i 6 aZ - 



(5) 
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For small e, i.e., close to the critical threshold distribu- 
tion, we expand 



ez + ln(l - ez) = -\e 2 z 2 - \t 3 z 3 
with the result 

D{A) _ A A - 1 e- A e 2 
N ~ 2A! 

r4 _A C 2 2, 



(9) 
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By use of Stirling approximation 

A! ~ A A e~ A V2^A , 



(10) 



(11) 



— a reasonable approximation even for small A — this 
may be written 
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We see from (11211 that there is a crossover at a burst 
length around A c , so that 



D(A) f (S/n) 1 /^ 2 A" 3 / 2 for A « A c 
N ~{ (8tt -V2a- 5 /2 forA>A c . 



(14) 
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FIG. 3: Distribution of first burst for the critical strength dis- 
tribution (Eqn.[5J with xo = 0.50a; m . The simulation results 
(squares) are based on 10 6 samples with N — 50000 fibers. 
The 'star' symbol stands for the analytic result (Ean. IT%jl . 



For xq — 0.45x m , we have A c = 50, so the final asymp- 
totic behavior is not visible in Fig. ^ The crossover is 
seen better for xq — 0.40x m . In Fig. [21 there is clearly a 
crossover near A = A c = 12.5. 

The phenomenon is not limited to the uniform thresh- 
old distribution. The £ = 3/2 power law 5] in the 
burst size distribution will appear whenever a thresh- 
old distribution is non-critical, but close to criticality. 
This can be seen from the expression for the average 
force F(x) oc xQ(x) on the bundle. The critical value 
x c corresponds to the maximum of F(x), which gives 
Q(x c ) = x c p{x c ). For the weakest fiber strength Xo very 
close to x c and A finite, the integrand in (4) therefore 
approach a constant times e~ A (x c — x). Then D(A) is 
proportional to A A ~ 1 e~ A /A! ~ (27r)" 1 / 2 A" 3 / 2 . On the 
other hand, when A ^> (x c — xq)~ 2 the generic asymp- 
totic behavior D oc A~ 5 / 2 will follow. 

Precisely at criticality A c = oo and the £ = 5/2 power 
law is no longer present. We can demonstrate by a ran- 
dom walk argument that at criticality the burst distri- 
bution follows a 3/2 power law. The load on the bundle 
when the fcth fiber with strength Xk is about to fail is 



F k = Q(x k )x k , 



(15) 



and the fluctuations of this load determine the size of 
the bursts. The probability p(f) df that the difference 
Fk+\ — Fk is in the interval (/, / + df) has been shown 
to be 
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gK) e -(f+Xk)p(x k )/Q(x k ) for J > _ Xk 

for / < -x ka . 



At criticality xp{x) — Q(x), which gives 

-l e -i-f/x for j > _ x 
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for / < —x 



(16) 
(17) 



FIG. 4: The burst distribution in the fuse model. System 
size is 100 x 100 and the number of samples was 300. On 
the average, catastrophic failure sets in after 2097 fuses have 
blown. The circles denote the burst distribution measured 
throughout the entire breakdown process. The squares de- 
note the burst histrogram based on those appearing after the 
first 1000 fuses have blown. The triangles denote the burst 
histogram after 2090 fuses have blown. The two straight lines 
indicate power laws with exponents -3.0 and -2.0 respectively. 



This can be considered as the step probability in a ran- 
dom walk. The random walk is unsymmetrical, but it is 
unbiased, < f >= 0, as it should be at criticality. 

Using the step probability (|17f) this may be evaluated 
61, with the result 
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The simulation results in Fig. [3] are in excellent agree- 
ment with this distribution. At the completion of a burst 
the force, i.e., the excursion of the random walk, is larger 
than all previous values. Therefore one may use this 
point as a new starting point to find, by the same cal- 
culation, the distribution of the next burst, etc. Conse- 
quently the complete burst distribution at criticality is 
proportional to l(T%Jl . i.e. essentially oc A" 3 / 2 . 

In order to test this crossover phenomenon in a more 
complex situation, we have studied the burst distribu- 
tion in the fuse model 0. The fuse model consists of 
a lattice where each bond is an ohmic resistor as long 
as the electrical current it carries is below a threshold 
value. If the threshold is passed, the fuse burns out ir- 
reversibly. The threshold t of each bond is drawn from 
an uncorrelated distribution p(t). The lattice is placed 
between electrical bus bars and an increasing current is 
passed through it. Numerically, the Kirchhoff equations 
are solved with a voltage difference between the bus bars 
set to unity. The ratio between current ij and threshold 
tj for each bond j is calculated and the bond having the 
largest value, max j{ij/tj) is identified and subsequently 
irreversibly removed. The lattice is a two-dimensional 
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square one placed at 45° with regards to the bus bars. 
The threshold distribution is uniform on the unit interval. 
All fuses have the same resistance. The burst distribution 
follows the power law with £ — 3.0, which is consis- 
tent with the value determined by Hansen and Hcmmcr 
@ . We show the histogram in Fig. With a system 
size of 100 x 100, 2097 fuses blow on the average before 
catastrophic failure sets in. When measuring the burst 
distribution only after the first 2090 fuses have blown, a 
different power law is found, this time with £ = 2.0. After 
1000 blown fuses, on the other hand, £ remains the same 
as for the histogram recording the entire failure process, 
see Fig. 0] 



In conclusion, we have studied the burst distribution 
in the fiber bundle model and shown that close to catas- 
trophic failure it exibits a crossover behavior with two 
power law with exponents -3/2 and -5/2. In the critical 
situation a random walk argument gives a single power 
law with exponents -3/2. We show numerically, that 
the same crossover — but with different values for the 
exponents — in the two-dimensional fuse model. This 
crossover signals that catastrophic failure is imminent. 
This signal does not hinge on observing rare events, and 
therefore has a strong potential as a useful detection tool. 
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